Dipole resonances in light neutron-rich nuclei studied with time-dependent 
calculations of antisymmetrized molecular dynamics 
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In order to study isovector dipole response of neutron-rich nuclei, we have applied a time- 
dependent method of antisymmetrized molecular dynamics. The dipole resonances in Be, B and C 
isotopes have been investigated. In 10 Be, 15 B, 16 C, collective modes of the vibration between a core 
and valence neutrons cause soft resonances at the excitation energy E x = 10 — 15 MeV below the 
giant dipole resonance(GDR). In 16 C, we found that a remarkable peak at E x — 14 MeV corresponds 
to coherent motion of four valence neutrons against a 12 C core, while the GDR arises from the core 
vibration in the E x > 20 MeV region. In 17 B and 18 C, the dipole strengths in the low energy region 
decline compared with those in 15 B and 16 C. We also discuss the energy weighted sum rule for the 
El transitions. 



I. INTRODUCTION 



in 
o 
o 

w>: 

In neutron-rich nuclei, there often appear exotic phenomena which are very different from those in stable nuclei 
due to the excess neutrons. Some properties of neutron-rich nuclei are concerned with differences between proton and 
neutron densities. Neutron halo and skin structures are typical examples. Another related subject is deformations of 
proton and neutron densities. For example, the difference of the deformations between proton and neutron densities 
were theoretically suggested in some Be, B and C isotopes [1-3]. These phenomena imply that the structure of the 
f* — ' nuclei far from the /3-stability line often contradicts to the usual understanding for stable nuclei where the proton and 
in > neutron densities are consistent with each other in a nucleus. They may lead also to exotic phenomena in excitations 
and reactions. One of the current issues is the dipole excitations in neutron-rich nuclei [4-21]. Since it is quite 
natural to expect the low-energy isovector-dipole excitations due to the difference between the proton and neutron 
densities, interests are attracted to the soft resonances below the giant dipole resonance(GDR), their collectivity and 
the contributions of the excess neutrons. In fact, the features of dipole transitions in neutron-rich O isotopes have 
been studied experimentally [20] and theoretically [11,13,16]. They have been found to be different from those in the 
| stable nucleus, 16 0, especially in the low energy region below the GDR. The dipole excitations have been studied 
also in C isotopes by Suzuki et al. with shell model calculations [17], where they suggested that coherent neutron 
Q transitions enhance the strengths at the excitation energy E x = 10 — 15 MeV. 

Our present interest is in the isovector dipole excitations in the light neutron-rich nuclei and in the effect of 
the ground state properties like the deformations on the El strengths. A method of antisymmetrized molecular 
J> | dynamics(AMD) [1] is one of powerful approaches for nuclear structure study. The method is superior especially in 
the description of cluster aspect, which is important in light unstable nuclei as well as in light stable nuclei [22]. In 
the systematic studies of Be, B, and C isotopes performed with the AMD method, a variety of structure such as the 
neutron skin and the deformed states have been suggested in those nuclei, and some of them have been discussed 
in relation to cluster aspect [1,22]. The experimental data for various properties of the neutron- rich Be, B and C 
isotopes have been successfully reproduced by the AMD calculations. We should stress that the AMD calculations 
well agree to the experimental data of quadrupole moments and E2 strengths in neutron-rich B and C, which can 
not be reproduced by the shell model calculations without using system-dependent effective charges. For the study of 
dipole excitations in the AMD framework, we apply a time-dependent method and calculate the dipole strengths in a 
similar way to the the time-dependent Hartree-Fock(TDHF). The point is that we are able to study dipole resonances 
with the framework which can describe cluster aspect. One of the advantages of the time-dependent AMD is that we 
can link the excitations with such collective modes as core vibration, core-neutron motion and inter-cluster motion 
which should be important to understand the role of the excess neutrons in the dipole resonances. 

The time-dependent method of AMD have been proposed and applied to heavy-ion reactions by Ono et al. in 1992 
[23,24] earlier than the application of the AMD to the nuclear structure study. For collective motion on the static 
solution, however, the time-dependent AMD calculations have not been performed yet. In this paper, we formulate 
a method based of the time-dependent AMD for the study of El response in analogy to the TDHF. In order to see 
its validity, we first apply it to 12 C and 18 O, and show comparison with the experimental data and other theoretical 
calculations. Then we apply this method to Be, B and C isotopes and discuss the properties of dipole strengths in the 
neutron-rich nuclei. We try to see how the dipole strength distribution is influenced by such structure as deformations, 
neutron skin, the existence of core and clusters. 
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This paper is organized as follows. In the next section, we explain the formulation of the present method for the 
El response, which is based on the time-dependent AMD. Adopted effective nuclear forces are described in III. In 
IV, we show the results of 12 C, 18 and the comparison with the experimental data, and also the dipole transitions in 
Be, B, C isotopes. The discussion of the El excitations in neutron-rich Be, B, and C isotopes are given in V. Finally, 
in VI we give a summary. 



II. FORMULATION 



We explain the formulation of the time-dependent version of AMD for study of isovector dipole excitations. By 
simulating time evolution of the collective motion on the static solution with the time-dependent AMD, we can 
calculate the response of a nucleus to external dipole fields and obtain the dipole strengths in the similar way to 
TDHF approaches. 

The time-dependent method in the AMD framework is described in Refs. [23,24], where the method has been 
applied to heavy-ion collisions. Concerning nuclear structure study with AMD methods, the static version of AMD 
and its extended versions are reviewed in Refs. [22,25]. 



A. Wave function 



The wave function for a A-nucleon system(A is a mass number) is given by a single Slater determinant of Gaussian 
wave packets as, 

$ = -7=Mfl,f2,--,<fiA}, (1) 

where the z-th single-particle wave function is written as follows, 

<Pi = <PziXin, (2) 
<Az i (r J ) = (|) 3 / 4 exp{-Kr,--?|) 2 }, (3) 

Xi = (^+6)XT + (^-6)Xi- ( 4 ) 

Here, the spatial part of the i-the single particle wave function is given by a located Gaussian wave packet, whose 
center is represented by the complex parameter, Zj. The parameter £j indicates the orientation of the intrinsic spin, 
and the iso-spin function is up(proton) or down(neutron). 

In the present work, the orientation of the intrinsic spin is fixed to be up or down as & = {1/2, —1/2} for simplicity. 
In this AMD wave function, all the centers of Gaussian wave packets for A nucleons are independent variational 
parameters, and a set of parameters Z = {Zi, Z 2 , • • ■, Za} specifies the total wave function <&(Z) of the state. This 
is the simplest version of AMD wave function, and the parity and angular momentum projections are not performed 
in the present work. 



B. Equation of motion 



In the time-dependent version of AMD, components of Z are considered to be time-dependent parameters as 
explained in [24]. The time evolution of the system is determined by the time-dependent variational principle, 
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This leads to the equation of motion with respect to Z, 

= dZ* 



ih^CjgjrZjr = ^ r and c.c, (6) 



3,r 

where a, r = x, y, z, and H is the expectation value of the Hamiltonian H, 
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(7) 



= ln($(Z)|$(Z)) (8) 

is a positive definite Hermitian matrix. These equations, (6), (7) and (8) are derived in general from the time- 
dependent variational principle for a given wave function parametrized by complex variational parameters. In case 
of the AMD framework, the time evolution of a system is described by the motion of the centers of Gaussian wave 
packets. 

Although stochastic collision process has been introduced in the study of heavy- ion collisions [24] , we do not put it 
in the present framework. 



C. Response to dipole fields 

In order to calculate the response to external fields, we first solve the static problem to obtain the optimum solution 
\E (0 for the ground state. We perform the energy variation of AMD wave function with respect to the variational 
parameter Z by using the frictional cooling method(a imaginary-time method) [22-24]. We obtain the optimum 
parameter Z°, which gives the energy minimum state "J" = $(Z°) in the AMD model space. Then, we boost the 
^° instantaneously at t=0 by imposing an external perturbative field, V ey ±{r,t) = eF(r)S(t), where e is an arbitrary 
small number. This results in an initial state of the time-dependent calculation as follows: 

*(t = 0+) = e- ttF ^° = e~ teF <!>(Z Q ) (9) 

In the calculation of El resonances, the external field is chosen to be the dipole field as, 

A 

F(r) = eM(El,n) = ^ e^r^r,), (10) 

i 

where e roc is the El recoil charge, Ne/A for protons and — Ze/A for neutrons. Then the initial state ^{t = 0+) is 
written with a single AMD wave function <f>(Z(t = 0+)) by simply transforming the parameters Z° = {Z°, Z°, • • • , Z A } 
as follows: 

Z i (t = O+) = Z?-f0ii, (11) 

where e M is the unit vector. Although an extra normalization factor of the wave function arises from this transfor- 
mation, it gives no effect on the physical quantities because the AMD framework is always based on the normalized 
wave functions. 

By using the equation of motion (Eq.(6)), we can calculate the time evolution of the system, = $(Z(t)), from 
the initial state *f>(t = 0+) = $(Z(t = 0+)) following the time-dependent AMD. The transition strength is obtained 
by Fourier transform of the expectation value of A4(E1, fi) as follows, 

dB(cj;El,n) _ ^ ..w™ Mnxi2c/ n 1 



1 f 00 

Y j \{n\M{El, l i)\0)\ 2 5{w-u n ) = Im / dt(V(t)\M(El, M )|*(i))e iwt , (12) 

„ 71-6 Jo 



where |0) is the ground state and |n) is the excited state with the excitation energy Uu) n . In the deformed nuclei, 
Eq. (12) gives the El transition strengths in the intrinsic state because the total angular momentum projection is 
not performed. Assuming the strong coupling scheme, we calculate the B{E1) in the laboratory frame by sum of the 
intrinsic El strengths as follows: 

dB(iu;El) _ ^ dB(u;El,K) 

In the practical calculation, we impose the El field with respect to each direction, x, y, z, independently, and sum 
up the strengths instead of the sum of K = 0, ±1. In the present framework, dB(u);El)/dw consists of discrete 
peaks in principle, because the present AMD is a bound state approximation and continuum states are not taken into 
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account. We introduce a smoothing parameter T, add an imaginary part iT/2 to the real excitation energy E x as 
Tluj = E. x + iT/2, and calculate the B{E1) with Eq.(12) by performing the integral up to finite time. This smoothing 
can be considered to simulate the escape and the spreading widths of the resonances. 
The photonuclear cross section a(u) is related to the transition strength B(u; El) as, 



, , 1671-3 dB(u;El) 
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(14) 



III. EFFECTIVE NUCLEAR INTERACTIONS 



We use an effective nuclear interaction which consists of the central force, the spin-orbit force and the Coulomb 
force. In the present work, we adopt MV1 force [26] as the central force. This force contains a zero-range three-body 
force in addition to the finite-range two-body interaction: 



yMVl = ^y(2) + V {3 \ 

i<j<k 



i<j 



= (w + bP a - hP T - mP a P T ) 
V (3) =t 3 5(r i -r j )8{r j -r k ), 



Vi exp 



r 2 \ f r 2 - 

13 1 + V 2 exp ' 13 
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at, 



(15) 

(16) 
(17) 



where P a and P T denote the spin and isospin exchange operators, respectively. The two-body part contains 
Wigncr(u) = 1 — m), Bartlett(fe), Hciscnbcrg(/i) and Majorana(m) terms. Concerning the spin-orbit force, the same 
form of the two-range Gaussian as the G3RS force [27] is adopted. Coulomb force is approximated by the sum of 
seven Gaussians. 

In the present work, we use the same interaction parameters as used in Ref. [1] except for 18 0. Namely, we use the 
case 3 of MV1 force and choose the Bartlett, Heisenberg and Majorana parameters as b = h = and to — 0.576, 
respectively. The strengths of the spin-orbit force is chosen as uj — —Uu = ui s = 900 MeV. For 18 0, we can not 
obtain a stable solution of the AMD wave function without parity projection in case of the parameter m = 0.576 due 
to a problem of the numerical calculation. It is because the Gaussian centers Zj gather to the origin and the norm of 
the AMD wave function becomes almost zero in the energy variation. In order to avoid this problem, we use a slightly 
large Majorana parameter as to — 0.62 instead of to = 0.576. We note that the properties of the ground state are not 
qualitatively unchanged in the parameter range to = 0.576 ~ 0.63 in most nuclei [22]. 



IV. RESULTS 

We apply the present method of AMD to dipole excitations in 8 > 10 > 14 Be, n.iS.^B, 12,16,18,20 C) 18q In thc present 
AMD method without parity and spin projections, we can not obtain static solutions for the N = 8 isotones by the 
cooling method due to the divergence of the inverse norm of the wave function, because a system with N — 8 favors 
the p-shell closed state, which is written by the AMD wave function with the zero-limit of Gaussian centers (Z) for 
all the neutrons. 



A. Properties of ground states 

The wave functions of the ground states^ ) are obtained by the energy variation for the AMD wave function 
without spin-parity projections. The width parameter v is fixed and chosen to be an optimum value for each nucleus 
to give the minimum energy of the ground state in most cases. For 15 B, 16 C, 18 C, and 18 0, we use a slightly larger 
width parameter than thc optimum value to avoid the numerical problem in thc norm of the wave function. The 
adopted v values are listed in Tabled 

As is suggested in Refs. [1-3,22], the shape of the proton and neutron density distribution rapidly changes with 
the variation of the proton and neutron numbers. The root-mean-square radii for the ground states^ ) are shown in 
Fig.l. In each series of isotopes, the neutron radius is enhanced in the neutron-rich nuclei with the increase of neutron 
number. In Fig. 2, we show the deformation parameters (/?, 7) for proton and neutron densities. The results are the 
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TABLE I. The adopted width parameters(^) of the AMD wave functions. 
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FIG. 1. The root-mean-square radii of the ground states obtained by the AMD calculation. The radii for the proton(neutron) 
density distributions of the ground state, are plotted by crosses(squares). The circles indicate the nuclear matter radii. 



qualitatively same as the previous results obtained by the simple AMD calculations [2,3]. The difference of the shape 
between protons and neutrons is found in the 7 parameter as well as the (3 value of some nuclei. The discrepancy of 7 
is remarkable in 10 Be and 16 C. Namely, opposite deformation between proton and neutron densities appears in these 
nuclei as discussed in Ref. [2,3]. 



B. Energy weighted sum rule 

The energy weighted sum rulc(EWSR) for isovector dipole resonances is given by 

s(Ei) = jn J B{ 2 E1) ^. (18) 

If the interaction commutes with the El operator, S(E1) is identical to the classical Thomas-Reiche-Kuhn(TRK) 
sum rule: 

9e 2 NZ 

where M is the nucleon mass. Due to the contributions of exchange terms and momentum dependent terms, the 
interaction is usually incommutable with the El operator and S(E1) is enhanced compared with S(TRK). Following 
the method explained in Ref. [28], we can calculate the EWSR with the initial-state expectation value of the double 
commutator of the Hamiltonian with the dipole operator F 

S(E1) = ^°\[F,[H,F]]\*°). (20) 

We estimate the enhancement factor, n = S(E1)/S(TRK) — 1, for the present interaction by neglecting the contribu- 
tion of the spin-orbit force. The incommutable terms in the present interaction come from Heisenberg and Majorana 
exchange terms in the two-body central force V^ 2 \ If we write the two-body force as = v(rij) + v T (rij)Ti ■ Tj (r 
is the isospin SU(2) generator), the enhancement AS(E1) = S(E1) — S(TRK) is given as follows [28], 

^S(El) = -l-e 2 ^ \J2r^(r ij ) (t x (i)t x (j) + t y (i)t y (j)) |*°). (21) 

Here t = \t. By calculating the expectation value, Eq.(21), for the static solution >F = ^(Z ), we can obtain the 
values S(E1) and k. 

In table II, the calculated S(E1) and k are shown. In order to demonstrate that the sum rule is kept in the present 
framework of the time-dependent AMD, we compare the value S(E1) = S(TRK) + AS(E1) given by the static 
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FIG. 2. The deformations of the ground states calculated by AMD. The deformation parameters (a)/3 and (b)7 for the 
proton (neutron) density distributions in *° are plotted by crosses(squares). Circles are for the deformations of nuclear matter 
density. 



calculations with the EWSR value S(E1; total) obtained by integrating the strengths Eq.(18) calculated with the 
time-dependent AMD. As is shown in table II, S(E1) = S(E1; total) is practically satisfied. It is reasonable because 
the present calculation is regarded as a method based on the small-amplitude TDHF. The enhancement factor k is 0.71 
and 0.74 for 12 C and 18 0, respectively, and it is k = 0.6 ~ 0.7 for the neutron-rich Be, B, and C isotopes. These are 
consistent with k — 0.4 <~ 0.8 due to the effect of the exchange mixtures of two-body interactions stated in Ref. [29]. 
In the shell model calculations [13,17], the values k = 40% - 50% for 12 C and 18 are obtained for S(E1;E X < 40 
MeV) integrated up to 40 MeV, while k = 0.13 for S(E1; E x < 40 MeV) in 18 is obtained by quasiparticle random 
phase approximation(QRPA)+phononcoupling model [16]. In the experimental photonuclear reactions, the observed 
cross section integrated up to 30 MeV for 12 C [30] exhausts 63% of TRK sum rule value, and 90 % of S(TRK) is 
exhausted by EWSR integrated up to 42 MeV in 18 [31]. As shown later, the EWSR is dominated by the GDR in 
the present results, which means that the calculated GDR should be quenched and the large fraction of the strength 
should be in the higher energy region than the GDR. 



C. Dipole resonances 



1. 12 C and 18 O 



We first show the results of the dipole resonances in 12 C and 18 O, and compare the results with other theoretical 
calculations and experimental data to see validity of the present method. The photonuclear cross section of 12 C and 
18 is plotted as a function of the excitation energy in Fig. 3 and 4. Thin dash-dotted, solid, and dotted lines indicate 
the contribution of vibration for the x,y, and z-directions, respectively. Here and hereafter, we chose the x, y, and 
z axis as (^ \x 2 \^ ) < (*°|y 2 |*°) < ('5 |z 2 |* ) and (xy) = (yz) = (zx) = 0. The thick solid lines correspond to 
the total strengths. We use the smoothing parameter T = 1,2,4 MeV. In the present results, the GDR peak lies at 
E x = 26 MeV and E x = 28 MeV in 12 C and 18 0, respectively. These peak positions are about 4 MeV higher than the 
observed GDR peaks [30-33], and also higher than other theoretical values of the shell model [13,17] and the QRPA 
calculations [16]. Compared with the observed photonuclear cross section, we need a smoothing parameter T > 4 
MeV to reproduce the width of the GDR. The reason for such a large T is considered to be due to the limitation of 
the present model space and lack of the effects of continuum states. For the quantitative discussion of the magnitude 
of the GDR strength, further quenching and the spreading are needed in the present calculations. 

Although the quantitative description of the peak positions and the magnitudes of the GDR are not sufficient, 
the characteristic behavior of the calculated cross section is in reasonable agreement with that of the experimental 
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TABLE II. The energy weighted sum rule(EWSR) values of El transitions. The S(E1; total) values are obtained by in- 
tegrated the strengths up to a enough large energy E x = 100 MeV in the time-dependent calculations of the AMD. The 
enhancement AS(E1) — S(E1) — S(TRK) is given by the ground state expectation values of the double commutator, Eq.(21). 
Here we ignore the contribution of the spin-orbit force. The values of the enhancement factor k = S(E1)/S(TRK) — 1 are also 
shown. The unit is MeV e 2 fm 2 for the energy weighted sum rule values. 





z>(£jl; total) 


O { 1 KA J 


Zlo(iil ) 


S(E1) 


i>\EL, total] / b [hjL ) 

w 

70 


K 


8 Be 


51 


30 


21 


51 


100 


0.71 


10 Be 


58 


36 


23 


58 


99 


0.63 


14 Be 


69 


42 


27 


69 


99 


0.62 


n B 


67 


40 


27 


67 


99 


0.66 


15 B 


81 


49 


32 


81 


100 


0.63 


17 B 


85 


52 


33 


86 


100 


0.63 


12 C 


76 


44 


32 


76 


100 


0.71 


ie c 


92 


56 


36 


92 


100 


0.65 


is c 


98 


59 


39 


98 


100 


0.65 


20 C 


103 


62 


40 


102 


101 


0.66 


18 Q 


115 


66 


49 


115 


100 


0.74 



data, and is consistent with other theoretical calculations. Since C has the oblate deformed ground state, the 
vibration for the y and z-axes forms the GDR in the same energy, which results in an enhancement of the lower part 
of the GDR. In the results of 18 O, significant dipole strengths are distributed in the energy region below the GDR 
due to the valence neutrons. The strong resonances have been experimentally observed in the region 10 — 15 MeV 
[20,31], and about 8% [20] of the TRK sum rule is exhausted by the integrated strength of the experimental data 
up to E x = 15 MeV. These low-lying resonances arc well described by the shell model calculations [13], which gives 
S(E1; E x < lbMcV) /S(TRK) = 6%. In the present results, the excitation energies of these low-lying resonances seem 
to be overestimated compared with those of the shell model and QRPA calculations as well as the GDR. Namely, the 
strengths of the low-lying peaks are distributed in the E x = 15 ~ 20 region, and S(E1;E X < 17 MeV)/S(TRK) = 6% 
and S(E1;E X < 15MeV)/ S(TRK) = 3% in the present results. Considering the shift of the excitation energies, we 
can state that the calculated strengths of the low-lying resonances reasonably agree to the experimental data. 

One of the reason for the overestimating excitation energies of the dipole resonances is considered to be because 
the surface diffuseness may be underestimated by the simple AMD wave function. It might be improved by the 
extension of the model wave function such as the deformed base AMD proposed by one(M.K.) of the authors and his 
collaborators [34]. 

2. C, B and Be isotopes 

Next we investigate the dipole resonances in neutron-rich Be, B, and C isotopes. We show the El strengths in 
8 Be, 10 Bc, and 14 Be in Fig. 5, and the photonuclear cross section in the left column of Fig. 8. In 10 Be, the dipole 
resonances in 10 Be can be decomposed into two parts E x < 15 MeV and E x > 20 MeV. The former consists of the 
soft resonances with the dominant strengths in 10 < E x < 15 MeV. The latter contains the GDR with a double peak 
structure with 7 — 8 MeV energy splitting, which is similar to that of 8 Be. In the E x < 5 MeV region, we find a peak 
with the strength B(E1) = 0.06 e 2 fm 2 . We consider this is a 1~ state and corresponds to the known 1~ state at 5.96 
MeV. The present low-lying peaks in the 10 < E x < 15 MeV originate in the cluster structure. The details will be 
discussed in the next section. In the TDHF+absorbing boundary condition(ABC) calculations, there are not such the 
significant strengths of the soft El resonances [18]. On the other hand, the GDR of the TDHF+ABC calculations 
is consisitent with the present results. In 14 Bc, the GDR splits into two peaks at E x = 15 MeV and at E x > 25 
MeV. The lower peak appears in the vibration along the longitudinal axis(z). As seen in Figs.l and 2, 14 Be has the 
large prolate deformation of the neutron density as well as the large neutron radius, and hence, it has the enhanced 
neutron skin structure along the longitudinal direction. The decrease of the excitation energy of the lower GDR peak 
is naturally understood because of the developed neutron skin. Also in the TDHF+ABC calculations [18], the GDR 
for the longitudinal motion appears at E x = 15 MeV, while the higher peak for the transverse motion is around 
E x = 25 MeV. Although the peak position of the GDR for each direction is similar to the present results, the GDR is 
not splitting in the TDHF+ABC results because the widths are largely spread. Another difference with the present 
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FIG. 3. The calculated photonuclear cross section of C. The smoothing parameter F = 1, 2, 4 MeV are used. Thin 
dash-dotted, solid, and dotted lines are the contribution of vibration for the x,y, and z-directions, respectively. The thick solid 
lines indicate the total strengths. 
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FIG. 4. The calculated photo-nuclear cross section of O. The smoothing parameter F = 1, 2, 4 MeV are used. Thin 
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results is that there exists a very soft resonance at E x <~ 5 MeV in TDHF+ABC. These differences seem natural 
because of the following reason. The spreading and quenching of the GDR may be large in 1 Be, which has a small 
neutron separation energy and therefore the dipole strengths should be affected by the continuum states and the long 
tail of the neutron halo structure. These effects are not taken into account in the present framework, while they are 
included in the TDHF+ABC. In the shell model calculations, the large El strength is found in the low energy region 
(5 < E x < 15 MeV) when including the fp shell configurations with a Warburton-Brown(WBP) interaction [15]. 

The El strengths and the photonuclear cross section in B isotopes are shown in Figs. 6 and 8. In the B isotopes, 
the feature of the GDR changes reflecting the variation of the deformation as the neutron number increases. In 15 B, 
two peaks of the GDR appear at E x = 20 and 27 MeV (see Fig.6(b)). Due to the prolate deformation, the lower 
GDR at E x = 20 for the longitudinal(z-direction) vibration has a smaller strength than that of the higher one. The 
excitation energy of the GDR is smaller than that of the n B. One of the unique features in 15 B is that a soft resonance 
appears in E x = 10 — 17 MeV region, which exhausts 16% of the S(TRK) value. This soft resonance arises from 
the longitudinal vibration and decouples energ etically with the GDR region^ > 17) MeV. In 17 B, the GDR peaks 
spread over a wide energy region due to the triaxial deformation. The peak position of the lowest GDR further shifts 
toward the low energy region: E x ~ 18 MeV. We can not find strong soft El resonance other than the GDR in 17 B. 

We show the results of C isotopes in Figs. 7 and 8. In comparison between B(Fig.6) and C isotopes(Fig.7), it is 
found that the feature of the dipole transitions in 16 C is quite similar to that in 15 B, which has the same neutron 
number(7V = 10) with 16 C. Namely, the dipole strength for the longitudinals-axis) vibration splits into two peaks, 
the GDR at E x = 22 MeV and a soft resonance at E x = 14 MeV. As a result, 16 C has a significant dipole strength in 
the low-energy region(i? 2: < 17 MeV) below the GDR region. This soft dipole resonance at E x = 14 MeV is consistent 
with the shell model calculations [17], where a remarkable peak is found at the same energy. In 18 C, since it has a 
triaxial deformation as well as 17 B, the shape of the strength function in the GDR region is similar to that of 17 B, 
though the peak positions are slightly higher than those of 17 B. A difference between 18 C and 17 B is the soft dipole 
strengths in the energy region E x < 17 MeV. Although there is no noticeable peak in this energy region, we find 
some fractions of the dipole strengths in 18 C rather than in 17 B. In 20 C with an oblate deformation, the shape of the 
strength function dB{El,uu) / duj is similar to the 12 C, while the peak positions are 4-5 MeV lower than those of 12 C. 



V. DISCUSSION 



As shown before, the remarkable peaks of the soft resonances are found in the dipole strengths of 10 Be, 15 B and 
16 C. They are separated clearly from the GDR region. It is natural to expect that these soft resonances arise from 
the coherent excitations of excess neutrons. In order to link the El resonances with collective motions we analyze the 
time evolution of the single-particle wave functions. In the time-dependent AMD, the expectation value of the dipole 
operator M(E1, [i = 0) for the ^(t) = $(Z(t)) is directly related to the real part of the centers of the single-particle 
Gaussian wave packets 



(*(t)|M(£l,/i = 0)|*(t)) 



N 



z iz {t) 



z 

A 



i=Z+l 



z iz {t) 



(22) 



where the i-th particle is a proton(neutron) for i = 1, • • • , Z(Z + 1, ■ ■ ■ , A), and Z iz is the z-component of the center Z^ 
for the i-th single particle Gaussian wave function. It should be stressed that the El excitations are expressed by the 
motion of the centers of single-particle Gaussian wave packets. Since the El strength is given by the Fourier transform 
of Eq.(22) as explained in (II), we can examine separately the contribution of the motion of each single-particle wave 
packet to the dipole strengths by Fourier transform of Re [Z iz (t)/ ' \fv\ and explicitly see collective modes. As discussed 
later, in case that a collective mode due to the inter-cluster motion appears, the mode can be seen as a peak at the 
corresponding excitation energy in the sum of the components for nuclcons in each cluster, 



— Im ( dtY Rc[Z lz /^]e luJt , 



(23) 



where C\, C2, • • • are the constituent clusters, and e is the same parameter in Eq.(12). 

In Fig. 9, we illustrate the density distribution and the spatial configuration of the Gaussian centers Re[Zj/-^/77] in 
the static solution of 16 C. There is the difference between proton and neutron densities in the ground state. The El 
transitions is described by the small-amplitude motion around this static solution. We see a 2n+ 12 C+2n configuration 
in the spacial configuration of the Gaussian centers, which forms the prolate neutron deformation with the longitudinal 
z axis. After the instantaneous external dipole field M.(E\,n = 0) is imposed, four valence neutrons coherently move 
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FIG. 5. The calculated El transition strengths of 8 Be, 10 Be, and 14 Be. The smoothing parameter is chosen to be F — 2 
MeV. Thin dash-dotted, solid, and dotted lines are the contribution of vibration for the x,y, and z-directions, respectively. The 
total strengths are shown by the thick solid lines. 
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FIG. 6. The calculated El transition strengths of n B, 15 B, and 17 B. The smoothing parameter is chosen to be F = 2 MeV. 
Thin dash-dotted, solid, and dotted lines are the contribution of vibration for the x,y, and z-directions, respectively. The total 
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against the core 12 C with the oscillation energy E x <~ 14 MeV to form the soft dipole peak. On the other hand, it is 
found that the strengths in the GDR region (E x > 20 MeV) arise from the motion of the nucleons within the 12 C core. 
In figure 11(a), we show the strengths of the motion for 4 valence neutrons, 6 protons and 6 neutrons in the 12 C core. 
It is found that the valence neutrons move with a negative strengths against 6 protons and 6 neutrons in the E x ~ 14 
MeV region, while, in the GDR region, the dipole strength is dominated by the relative motion between 6 protons 
and 6 neutrons inside the 12 C core. The reason why four neutrons move coherently in the 2n+ 12 C+2n configuration 
is easily understood as follows. Since the 2n+ 12 C+2n configuration is linear, one can consider a configuration with 
two dineutrons on the opposite sides of the core( 12 C). Let us imagine the inert three clusters(2n+ 12 C+2n), which 
are connected with two identical springs as shown in Fig. 9(d). In the motion along the longitudinal axis, there are 
two eigen modes of the oscillation. One is the mode where two dineutrons move in phase and the core moves in 
the opposite way, and the other is the one with the opposite motion of the dineutrons to each other. The former 
corresponds to the isovector dipole mode. Thus, the coherent motion of the valence neutron can be interpreted by the 
relative motion in the 2n+ 12 C+2ra configuration. The soft dipole resonance due to the excess neutrons in the present 
result of 16 C well corresponds to the shell model calculations [13], where the Op — > ls0<i and IsOd —* 0/lp transitions 
work coherently to enhance the strength at E x — 12 — 14 McV in 16 C. 

We show the density distribution of the ground state of 10 Be in Fig. 10. In order to understand collective motion 
for the soft dipole resonances, it is useful to regard the 10 Be as the a+ 6 He cluster state. In the analysis of motion of 
the Gaussian centers, it is found that the strength at E x = 10—15 MeV contains two independent modes. One is 
the inter-cluster motion between a and 6 He, which contributes to the resonance at E x = 12 MeV in the longitudinal 
vibration along the z-axis. The other is the coherent motion of the valence neutron against the core 8 Be, which 
results in the resonance at E x = 14 MeV in the vibration along the y-axis. In Fig. 11(c), we show the strength of 
the a- 6 He inter-cluster motion. It has a dominant peak at E x = 10 — 15 MeV, which corresponds to the soft peak 
in the z-component. On the other hand, the GDR is described by the motion inside the core 8 Be. Also in 15 B, we 
find that the coherent motion of the valence neutrons contributes to enhance the strengths of the soft resonance. It 
is concluded that the remarkable peaks at E x = 10 — 15 MeV in 16 C, 10 Be and 15 B arise from the coherent motion of 
the valence nucleons, which decouples with the motion inside the core. 

We show the calculated photonuclear cross section in Fig. 8. The shape of the strength function in the GDR 
region^a; > 17 MeV) has a close relation with the deformation of the system. In the oblately deformed system such 
as n B and 20 C as well as 12 C, the GDR splits in two part. The lower GDR peak has large transition strength. In the 
neutron-rich nucleus, 20 C, the excitation energy of the GDR is the lowest among the three nuclei. Also in the 8 Be and 
14 Be with the prolate deformation, the GDR splits into two, but the higher GDR has larger strength in contrast to 
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FIG. 9. The density distribution of the ground state of 16 C. The densities are integrated along the z-axis. (a)Matter density, 
(b) proton density, (c) neutron density are shown. The spatial configuration of the centers of single-particle Gaussian wave 
packets for protons and neutrons are also plotted by open circles in (b) and (c), respectively. Each circle consists of a spin- up 
nucleon and a spin-down nucleon. (d): The schematic figure for the configurations of the Gaussian centers in 16 C written by 
2n+ 12 C+2n. 
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FIG. 10. The density distribution of the ground state of 10 Be. The densities are integrated along the x-axis. (a)Matter 
density, (b) proton density, (c) neutron density are shown. The spatial configuration of the centers of single-particle Gaussian 
wave packets for protons and neutrons are also plotted by open circles in (b) and (c), respectively. Each circle consists of a 
spin-up nucleon and a spin-down nucleon. (d): The schematic figure for the structure of Be with a+ He configuration. 
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FIG. 11. Strength of the motion of the Gaussian centers in the response of 16 C and 10 Be with an instantaneous dipole field 
along z-axis. The strength is given by Fourier transform of positions of the Gaussian centers as shown in Eq.(23). We divide 
A nucleons into some groups: C\, C2, ■ ■ ■ and sum up the components for nucleons in each group. (a)Motion of the Gaussian 



centers in 16 C. Ci consists of four valence neutrons, and (72(63) contains 6 protons(neutrons) in the 12 C core. (b)Motion in 
10 Be. Ci, C2 and C3 consist of two valence neutrons, 4 protons, and 4 residual neutrons. (c)Relative motion between 4 He(Ci) 
and 6 He(C2) clusters in 10 Be. The dotted, dash-dotted, and solid lines correspond to Ci, C2, and C3. 



16 



the oblate deformation. The lower resonance in 14 Be exists at E x = 15 McV. In triaxial nuclei such as 17 B and 18 C, 
the GDR consists of three peaks with 5 MeV of the energy splitting. However, considering spreading of the widths, 
three peaks may overlap to form a broad structure in the E x = 15 — 30 MeV. 10 Be and 16 C have the ground state with 
opposite deformations between proton and neutron densities. In spite of this unusual properties of the deformations, 
we can not find an abnormal feature in the dipole strength of the GDR region in these nuclei. As explained before, 
the remarkable soft peaks appear at E x = 10—15 McV due to the motion of the valence nucleons against to the core, 
while the GDR at E x > 17 McV arises from the vibration within the core. Therefore it is considered that the GDR 
may reflect mainly the features of the core nuclei instead of the deformations of the total system. In fact, the GDR 
of 10 Be lies at a similar excitation energy to that of 8 Be. 

Below the GDR region, we find the significant strengths of the soft resonances in such neutron-rich nuclei as 10 Bc, 
15 B, 16 C and 18 C due to the excess neutrons. Especially, 10 Be, 15 B and 16 C have the remarkable soft peaks, which 
are decoupled with the GDR. In the neutron-rich nuclei, the cluster sum rule [35] is a convenient measure to estimate 
the contribution of the motion of the excess neutrons in the dipole strengths [36] . Assuming clustering with a core 
and N v valence neutrons, we consider the core cluster with Z\ = Z,N\ = N — N v and the valence cluster with 
Z 2 = 0, N 2 = N v . The cluster sum rule is given as, 

h 2 9 (Z,N 2 - Z 2 N,) 2 2 h 2 9 Z 2 N 2 
clust 2m 4tt A(Z x + N 1 )(Z 2 + N 2 f 2m At: A(A - N V )N V 6 1 ' 

This value is the remainder when one subtract the core contribution of the classical EWSR from the total S(TRK) 
value. Consequently, S c \ us t is the margin which indicates the contribution of the excess neutrons. The integrated 
strength of the low-energy resonances should be compared with S c \ us t to see the softness and collectivity of the 
resonances due to the valence neutron motion against the core. In table III, the EWSR for the low-lying resonances 
are listed with the values of the classical TRK sum rule and the cluster sum rule (Schist)- I n the derivation of the cluster 
sum rule Schist, the core cluster is assumed to be 8 Be, 11 B, 12 C and 16 0, in Be, B, C and O isotopes, respectively. We 
show the EWSR values integrated up to E x < 15 MeV and E x < 17 MeV in the present results, and the EWSR with 
other calculations. The ratios of the EWSR to S(TRK) and 5 c i us t are shown in Figs. 12(a) and 12(b), respectively. In 
C isotopes, S(E1; E x < 17 MeV) is the largest in 16 C and it declines in further neutron-rich C isotopes. The present 
results of C isotopes well agree to the shell model calculations [13]. Also in B isotopes, the similar feature is found. 
Namely, S(E1;E X < 17 MeV) is the largest in 15 B and it decreases in further neutron-rich nucleus 17 B. The striking 
point is that the EWSR for the low-lying resonances is remarkably enhanced in the moderately neutron-rich nuclei with 
an appropriate number of excess neutrons, but it is suppressed in very neutron-rich nuclei. It is reasonable because 
the enhancement of the soft dipole strengths is due to the coherent motion of the valence neutrons relative to the 
core. It means that the decoupled collective modes appear based on the relative motion between the core and valence 
neutrons and the motion inside the core, in 15 B and 16 C. On the other hand, as the neutron skin develops in further 
neutron-rich nuclei N > 10, the motion of the excess neutrons is not decoupled but they join the neutrons inside the 
core. As a result, in 17 B and 20 C, the soft dipole mode is assimilated into the GDR, and the excitation energy of 
the GDR decreases. Also in 10 Be, the EWSR for the low-lying resonances is significant as well as 15 B and 16 C. In 
these nuclei, the cluster sum rule value 5 c i us t is almost exhausted by the calculated S(E1;E X < 17 MeV). In 14 Be, 
the S(El;E x < 17 MeV) is very large due to the peak at E x = 16 MeV. The reason for the enhanced S(E1; E x < 17 
MeV) in 14 Be is different from other nuclei( 10 Be, 15 B and 16 C). In case of 14 Be, the enhanced S(E1;E X < 17 MeV) 
does not originate in the soft resonance decoupled from the GDR, but the GDR for the longitudinal vibration itself 
contributes the EWSR for low-energy region because it becomes soft due to the prolate deformation with a developed 
neutron skin structure. 



VI. SUMMARY 



We applied a method of the time-dependent AMD to dipole transitions in light neutron-rich nuclei. We investigated 
the El resonances in Be, B, and C isotopes. It was found that the remarkable peaks appear in 10 Be, 15 B, and 16 C at 
E x = 10 — 15 MeV which almost exhaust the values of the cluster sum rule. These soft dipole resonances arise from 
the relative motion between excess neutrons and the core, which is decoupled with the motion inside the core. In other 
words, these soft resonances appear due to the excitation of the excess neutrons around the rather hard core. This 
nature of the neutron excitation and the inert core may have a link with such ground-band properties of the 1 C as the 
unusually small B(E2\l\ — > 0+) [37], which has been recently investigated in theoretical and experimental studies 
[3,38-40]. In further neutron-rich B and C isotopes with N > 10, the strength for the soft dipole resonances declines 
compared with that of 15 B and 16 C. It is considered to be because the motion of the excess neutrons assimilates into 
the neutron motion within the core. As a result, the excitation energies of the GDR decrease with the enhancement 
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TABLE III. The Energy-weighted sum rule for the soft dipole resonances. The energy weighted sum integrated up to E x < 17 
MeV and E x < 15 MeV are shown in the fourth and fifth columns. The core cluster in the derivation of the cluster sum rule 
Sdust is 8 Be, 11 B, 12 C and 16 0, in Be, B, C and isotopes, respectively. The smoothing parameter is chosen to be T — 2.0 
MeV. The unit is e 2 fm 2 MeV. 
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of the neutron-skin. It is striking that the strength for the soft dipole resonances does not necessarily increase with 
the increase of the excess neutrons. Instead, the feature of the soft resonances rapidly changes depending on the 
proton and neutron numbers of the system. The key of the soft dipole resonance with a remarkable strength is how 
the coherent motion of the valence neutrons is decoupled with the motion inside the core. 

The present method based on the time-dependent AMD is regarded as a kind of small-amplitude time-dependent 
Hartree-Fock calculations within the AMD model space. The point of the method is that we are able to study dipole 
resonances with the framework which can describe cluster aspect. In the AMD approach, the dipole excitations are 
expressed by the motion of single-particle Gaussian wave packets, because the expectation value of the dipole operator 
is related directly to the centers of the Gaussian wave packets. One of the advantages of the time-dependent AMD is 
that we can link the excitations with such collective modes as core vibration, core-neutron motion and inter-cluster 
motion, which should be important to understand the role of the excess neutrons in the dipole resonances. 

Recently, extended methods of time-dependent mean-field calculations have been proposed and applied to the dipole 
transitions in neutron-rich nuclei. In the TDHF+ABC approach, which have been applied to deformed neutron-rich 
nuclei by Nakatsukasa and Yabana [18], the effects of continuum states are taken into account. Another method is the 
time-dependent density-matrix theory which has been applied to 22 O by incorporating two-body correlations [7]. In 
the present work, the contributions of continuum states are omitted, and the detailed descriptions of wave functions 
and two-body correlations should be insufficient, as the model space is a simple AMD wave function written by a 
Slater determinant of Gaussian wave packets. Therefore, we put an artificial smoothing parameter to simulate the 
width of the dipole resonances, because it is difficult to describe the escape and the spreading widths of the resonances 
in the present framework. Further extensions of the model must be essential to give quantitative discussions of the 
excitation energies and the strengths of the dipole resonances in nuclei near the drip line. It should be necessary to 
solve the remaining problem of the soft resonances in halo nuclei [4-6,12,14,19,21]. 

We comment that the usual AMD wave functions applied to the nuclear structure study are the advanced ones with 
some extensions such as the parity and spin projections, deformed Gaussian base and generator coordinate method 
[22,34,41-43], though the present AMD wave function is the simplest one with no extension. A combination of the 
time-dependent method and the extended AMD wave functions should be needed to include higher correlations beyond 
the present calculations, and also to obtain better description of the ground state properties. Moreover, interaction 
dependence of the dipole transitions is a remaining problem. 
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